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We study the evolution of a stochastic hehcal magnetic field generated in the early 
Universe after the electroweak phase transition, using standard magnetohydrody- 
namics (MHD). We find how the coherence length ^, magnetic energy Eyi and 
magnetic helicity H evolve with time. We show that the self-similarity of the 
magnetic power spectrum alone implies that ^ ~ t^/'^. This in turn implies that 
magnetic helicity decays as ~ t"^", and that the magnetic energy decays as 
-Em ~ t"'^-^"^", where s is inversely proportional to the magnetic Reynolds num- 
ber Reyi. These laws improve on several previous estimates. 
SUSX-TH-03-002 NORDITA-2003-5 AP 



1. Introduction 

Magnetic fields are found everywliere in the Universe, from planets to 
galaxy clusters^'^. On the galactic scale and above, the strength is of order 
a few /iGauss, maintained by dynamo action, with a characteristic timescale 
of roughly a rotation period, lO'^ yr. A seed is required to start the dynamo, 
and a simple calculation^ based on the age of a typical galaxy shows that 
the seed field must have been about 10~^° Gauss. 
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There are many ideas for the origin of seed field. Most conservatively, a 
Biermann battery operated at the era of reionisation"^ . More speculatively, 
a field could have been generated in the early Universe^, which would have 
given rise to stochastic, homogeneous and isotropic magnetic and velocity 
fields, characterized by their power spectra and initial length scales. In 
order to know the strength and coherence scale of this field today, we must 
study freely decaying magnetohydrodynamic (MHD) turbulence^'^. This 
talk reports earlier work^ where we have developed a new framework for 
understanding decaying 3D MHD turbulence, in the case where the fields 
are close to being maximally helical, as in the mechanisms proposed earlier^. 

Of particular interest arc how the coherence length ^, the magnetic en- 
ergy Em and magnetic helicity H evolve with time. Their growth and decay 
are thought to be power laws, resulting from an inverse cascade^, in which 
power is transferred locally in fc-spacc from small to large scales. This is 
not a new problem: there arc models both with^"'^^'^^'^"' and without ^^'^^'"'^^ 
helicity, and numerical simulations have also been performed^^'^. We have 
now developed^'^ what we believe is the definitive answer for the evolution 
of helical magnetic fields between the electroweak phase transition and the 
time of e'^e~ annihilation. 



2. 3D MHD simulations of decaying turbulence 

The MHD equations in an expanding Universe are most conveniently ex- 
pressed in terms of conformally rescaled fields B (where B = V x A is 

the magnetic field in terms of the magnetic vector potential A), u (the 
velocity) and dissipation parameters v (kinematic viscosity), rj (magnetic 
diffusivity)^^'^. The current density J = V x B. 

In the ideal limit u = ij = 0, there is an important conserved quantity in 
addition to the total energy, which is the magnetic helicity H = / A - B d'^x, 
also known as the Abelian Chern-Simons number. Also important is the 
magnetic Reynolds number Reu = £,v/r], where ^ and v are the typical 
length scale and velocity of the system, because it measures the relative 
size of the non-linear term in the induction equation. 

We solve the MHD equations numerically^^, taking u and B to be ho- 
mogeneous and isotropic Gaussian random fields drawn from a power-law 
distribiition with a high wavcmmiber cut-off. The initial helicity, ranged 
from identically zero to maximal, where the inequality (fc)| < 2k^^Eyi[k) 
is saturated. Helicity is not exactly conserved in our simulations, and pro- 
viding that the magnetic energy spectrum Euik) decays faster than 
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we can define a helicity scale such that 

H = -2r,H/^l. (1) 

If wc assume that the evolution of is described by a power law ^ V" 
it is clear that only if r = 0.5 does the magnetic helicity show a power 
law decay H = Ho{t/to)~'^^ where s = (Cdiff/C-ff)^: in terms of the diffusion 
scale ^diff = 2TTy/rit. Additional length scales we consider are the integral 
scale = 27r J dkk^^ Eu{k) / J dkEM(k), the relative helicity scale = 
Tr\H\/EM and the magnetic Taylor microscale = 27ri?rms/'^rms) where 
firms and Jrms are the RMS magnetic field and current density respectively. 
It is plausible that all these scales are proportionally related and Fig. 1 
shows that this is indeed the case. From the bound on the helicity power 




Figure 1. Time evolution of the ratio of lengtli scale to the diffusion scale for one run. 
The notation is ^i/fdiff (dashed), ^H/Cdifl (dot-dashed), ^a/^difl (dotted) and ^T/Cdiff 
(continuous). The length scales are defined in Eq. (1) and the subsequent text. 

spectrum one can show J^rel = (^r/^i) < !> so if this bound remains 
approximately saturated, and the helicity scale goes as ^ t°'^, the decay 
law for the magnetic energy is 

i^M (2) 

Given Hrel = Cr/Cij it is seen from Fig. 1 that i?REL is indeed of order 
unity and does not decay markedly with time. 
To characterize the decay laws we define 

Q{t) = -tEu/Eu, R{t) = -tH/2H. (3) 

In Fig. 2a we have plotted R{t) versus the quantity s{t) = (^diff/Cff)^ 
for several runs with different initial conditions. This figtire tells us that 
the value of R is approximately independent of time which confirms the 
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Figure 2. (a, left) The (juaiitity R{t), defined by Eq. (3) versus the quantity s{t) = 
iidis/iH)^- The length scales and ^diff a^'e defined in Eq. (1) and the subsequent 
text, (b, right) Energy decay exponent Q, defined in Eq. (3), plotted against the inverse 
magnetic Reynolds number Re^ , calculated using the helicity scale ^h- The fit is a 
quadratic least squares fit using the final values. 

power law decay of H. Fig. 2a also indicates that the quantity s is also 
approximately independent of time, hence reinforcing the relation ~ 

Taking the relative helicity to be constant, it follows from the power-law 
behavior of H that Eq. (2) is the energy decay law. In the limit of exact 
conservation of magnetic helicity, s ^ 0, the magnetic energy must decay 
according to the power law Em ^ t~^'^. 

The parameter s has interesting physical significance. Note that if 
— vt, where v is the RMS velocity, (i.e. if the eddy turn-over time is t) 
then s ~ {2^)"^ / Rcm, where Rcm is the magnetic Reynolds number evalu- 
ated using the helicity scale ^h- We have measured the ratio / = vt/^u and 
Re^ for all runs, and find that they are both approximately constant, and 
are linearly related. This means there should be a linear relation between s 
and Re^, and hence a quadratic relation between Q{t), the energy decay 
exponent defined in Eq. (3), and Re^. Fig. 2b, showing Q and Re^, con- 
firms that this is indeed the case, with asymptote at large Rcm consistent 
with Q = 0.5. Finally, assuming only that the magnetic energy power spec- 
trum is self-similar, and that Ohmic dissipation, if not dominant, always 
contributes a constant fraction to the energy loss, one can show'' that the 
characteristic length scale of the field must scale as thus justifying 
the assumption made at the start of this section. 

3. Discussion and conclusions 

To summarize, we have studied the evolution of decaying 3D MHD tur- 
bulence involving maximally helical magnetic fields. For finite magnetic 
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difFusivity there emerges an important quantity s = (Cdiff /CH)^, where is 
the hehcity scale defined in Eq. (1), and ^diff is the diffusion scale. We find 
— vt, where v is the RMS velocity, and hence that s oc Re^, the mag- 
netic Reynolds number evaluated using the hclicity scale. The magnetic 
field coherence length (which can be equally well expressed as the inte- 
gral, helicity or relative helicity scales) goes as ^ ~ magnetic helicity 
Hm ^ t~^^ and magnetic energy Em ~ A corollary is that Reu is 

constant once the system has reached self-similarity. Furthermore, we can 
extrapolate to the limit of very large magnetic Reynolds numbers, useful 
for example in the early Universe, to find H constant and Em t~^-^. 
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